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ABSTRACT 

In earlier works we pointed out that the disk's surface layers are non-turbulent and thus highly con- 
ducting (or non-diffusive) because the hydrodynamic and/or magnetorotational (MRI) instabilities are 
suppressed high in the disk where the magnetic and radiation pressures are larger than the plasma 
thermal pressure. Here, we calculate the vertical profiles of the stationary accretion flows (with radial 
and azimuthal components), and the profiles of the large-scale, magnetic field taking into account the 
turbulent viscosity and diffusivity and the fact that the turbulence vanishes at the surface of the disk. 
Also, here we require that the radial accretion speed be zero at the disk's surface and we assume that the 
ratio of the turbulent viscosity to the turbulent magnetic diffusivity is of order unity. Thus at the disk's 
surface there are three boundary conditions. As a result, for a fixed dimensionless viscosity a-value, we 
find that there is a definite relation between the ratio 1Z of the accretion power going into magnetic disk 
winds to the viscous power dissipation and the midplane plasma-/3, which is the ratio of the plasma to 
magnetic pressure in the disk. 

For a specific disk model with 1Z of order unity we find that the critical value required for a stationary 
solution is f3 c s» 2Ar/(ah), where h the disk's half thickness. For weaker magnetic fields, (3 > f3 c , we 
argue that the poloidal field will advect outward while for (3 < /3 C it will advect inward. Alternatively, if 
the disk wind is negligible (TZ <C 1), there are stationary solutions with (3 ^> (3 C . 

Subject headings: accretion, accretion disks — galaxies: jets — magnetic fields — MHD — X-rays: 
binaries 



1. INTRODUCTION 

Analysis of the diffusion and advection of a large-scale 
magnetic field in a accretion disk with a turbulent viscos- 
ity and magnetic diffusivity arising from the magnetorota- 
tional instability (MRI) shows that a weak large-scale field 
diffuses outward rapidly (van Ballegooijen 1989; Lubow, 
Papaloizou, & Pringle 1994). We mention but do not con- 
sider here the opposite limit where the magnetic field is 
sufficiently strong that it suppresses the MRI instability 
so that the disk is non-turbulent but accretion occurs due 
to angular momentum outflow to a magnetic disk wind 
or jet (Lovelace, Romanova, & Newman 1994). Earlier, 
Bisnovatyi-Kogan and Lovelace (2007) pointed out that 
the disk's surface layers are highly conducting because the 
MRI instability is suppressed in this region where the mag- 
netic and radiative energy-densities are larger than the 
thermal gas energy-density. Rothstein and Lovelace (2008) 
analyzed this problem in further detail and discussed the 
connections with global and shearing box magnetohydro- 
dynamic (MHD) simulations of the MRI. Lovelace, Roth- 
stein, & Bisnovatyi-Kogan (2009; hereafter LRBK) devel- 
oped an analytic model for the vertical (z) profiles of the 
stationary accretion flows (with radial and azimuthal com- 
ponents), and the profiles of the large-scale, magnetic field 
taking into account the turbulent viscosity and diffusivity 
due to the MRI and the fact that the turbulence vanishes 
at the surface of the disk. 

Here, we require that the radial accretion speed be zero 
at the disk's surface, and we assume that the ratio of the 
turbulent viscosity to the turbulent magnetic diffusivity is 
of order unity as suggested by MHD shearing-box simula- 
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tions (Guan & Gammie 2009). For a fixed dimensionless 
viscosity a-value, we find that there is a definite relation 
between the ratio 1Z of the accretion power going into mag- 
netic disk winds to the viscous power dissipation and the 
midplane plasma-/?, which is the ratio of the plasma to 
magnetic pressure in the disk. 

Section 2 discusses the model for the flow and ordered 
magnetic field in a viscous diffusive disk. Section 3 dis- 
cusses the solutions for a specific disk model. Section 4 
gives the conclusions. 

2. THEORY 

Following LRBK we consider the non-ideal magneto- 
hydrodynamics of a thin axisymmetric, viscous, resistive 
disk threaded by a large-scale dipole-symmetry magnetic 
field B. We use a cylindrical (r, <f>, z) inertial coordi- 
nate system in which the time-averaged magnetic field is 
B = Brr+BcfjCb +B z z, and the time-averaged flow velocity 
is v = v r r + v<p4> + v z z. The main equations are 

= -Vp + pg + -3 x B + F v , 



dt 
dB 

~dt 



Vx(v x B) - V x (77V x B) 



(1) 

(2) 



These equations are supplemented by the continuity equa- 
tion, V ■ (pv) = 0, by V x B = 4ttJ/c, and by V • B = 0. 
Here, n is the magnetic diffusivity, F" = —V-T v is the vis- 
cous force with TJ k = —pi/(dvj/dxk+dvk/dxj — (2/3)5jkV- 
v) (in Cartesian coordinates), and v is the kinematic vis- 
cosity. For simplicity, in place of an energy equation we 
consider the adiabatic dependence p cx p 7 , with 7 the adi- 
abatic index. 
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We assume that both the viscosity and the diffusivity 
are due to magneto-rotational (MRI) turbulence in the 
disk so that 



a Tf 1 9(z) , 



(3) 



where V is the magnetic Prandtl number of the turbulence 
assumed a constant of order unity (Bisnovatyi-Kogan & 
Ruzmaikin 1976), a < 1 is the dimensionlcss Shakura- 
Sunyaev (1973) parameter, c s q is the midplane isothermal 
sound speed, fix = (GM/r 3 ) 1 / 2 is the Keplerian angu- 
lar velocity of the disk, and M is the mass of the central 
object. The function g(z) accounts for the absence of tur- 
bulence in the surface layer of the disk (Bisnovatyi-Kogan 
& Lovelace 2007; Rothstein & Lovelace 2008). In the body 
of the disk g = 1, whereas at the surface of the disk, at 
say zs, g tends over a short distance to a very small value 
<~ 10~ 8 , effectively zero, which is the ratio of the Spitzer 
diffusivity of the disk's surface layer to the turbulent dif- 
fusivity of the body of the disk. At the disk's surface the 
density is much smaller than its midplane value. 

We consider stationary solutions of equations (1) and (2) 
for a weak large-scale magnetic field. These can be greatly 
simplified for thin disks where the disk half-thickness, of 
the order of h = c s q/Q.k, is much less than r. Thus we 
have the small parameter 

e=^ = ^«l. (4) 

. r V K 

It is useful in the following to use the dimensionlcss height 
C = z/h. The midplane plasma beta is taken to be 

(5) 

where /3 = cl /v\ , vao = ^o/(47rp ) 1 ^ 2 is the midplane 
Alfven velocity. Note that the conventional definition of 
beta is 2/3. The rough condition for the MRI instability 
and the associated turbulence in the disk is /3 ^ 1 (Balbus 
& Hawley 1998) and this is assumed here. 

The three magnetic field components are assumed to 
be of comparable magnitude on the disk's surface, but 
B r = = on the midplane. On the other hand the 
axial magnetic field changes by only a small almount go- 
ing from the midplane to the surface, AB Z <~ eB r -C B z 
(from V • B = 0) so that B z w const inside the disk. As 
a consequence, the dBj/dr terms in the magnetic force 
in equation (1) can all be dropped in favor of the dBj/dz 
terms (with j = r, </>). It is important to keep in mind that 
Bj is the large scale field; the approximation does not ap- 
ply to the small-scale field which gives the viscosity and 
diffusivity. The three velocity components are assumed to 
satisfy v 2 z <C c 2 and v 2 <C v^. Consequently, v^,(r,z) is 

close in value to the Keplerian value Vx(r) = (GM/r) 1 / 2 , 
except in the outer disk layers where the radial magnetic 
force may be comparable with the centrifugal and grav- 
itational forces. We normalize the field components by 
B = B z (r,z = 0), with b r = B r /B , 6 = B^/ Bo, and 
b z = B z /B « 1. Also, we define = v<f,(r, z)/vx(r) and 
the accretion speed u r = — v r /(ac s o). For the assumed 
dipole field symmetry, b r and are odd functions of ( 
whereas u r and u§ are even functions. 

Integration over the vertical extent of the disk gives the 
average accretion speed 

2 &</,s+ 



u 



a/3£ 



(6) 



(LRBK) which is the sum of a viscous contribution, uq = 
3sk v (with fc„ a numerical constant of order unity) , and a 
magnetic contribution (cx b < ps+) due to the loss of angular 
momentum from the surface of the disk where necessarily 
b<j>s+ < (Lovelace, Romanova, & Newman 1994). 

We assume p cx p 7 so that the vertical hydrostatic equi- 
librium gives 



Po V 2 7 



1/(7-1) 



(7) 



for f) 3> 1. The density goes to zero at ( m — [27/(7 ~ 
l)] 1 / 2 . However, before this distance is reached the MRI 
turbulence is suppressed, and g(() in equation (3) is effec- 
tively zero. 

The different components of equations (1) and (2) can 
be combined (LRBK) to give the following equation for 
the radial accretion speed, 



,4«2 



d 2 



d 



.8/18 



- a 2 /3V 



d 2 



dc, 2 

- a 2 pV 



d 



d( 

d 2 (\d_ 
dC 2 \JdC 



du r 

_ d I u r 
oC, \pg 



du r 



u r 

pg 

+ zpv 2 ^- = 



(8) 



The equation can be integrated from ( = out to the 
surface of the disk £5 where boundary conditions apply. 

For specificity we take g(Q = (1 - CVCi) 5 , where 
Cs < Cm and S is a constant. That is, we neglect the ratio 
of the Spitzer diffusivity on the surface of the disk to its 
value in the central part of the disk. An estimate of Cs can 
be made by noting that /3(C) = 4irp(()/B 2 = /3(p) 7 w 1 at 
Cs. This gives Cl/Cl = l-/3- (7 ~ 1)/7 and Ps /p - /3" 1/7 - 

2.1. Boundary Conditions 

We consider only the solutions which have net mass ac- 
cretion, M = AirrhpoaCsoTi u r > 0, and have b^ < on 
the disk's surface. This condition on 6^5+ corresponds to 
an efflux of angular momentum and energy (or their ab- 
sence) from the disk to its corona rather than the reverse. 
The condition b c f,s+ < is the same as u r > uq, where 
m is the minimum (viscous) accretion speed. Note that 
1Z = u r /uo — 1 is the ratio of the accretion power which 
goes into the disk winds to the viscous power dissipation. 
Clearly, the condition on b c f,s+ implies that M > so that 
there is only one condition. In general there is a continuum 
of values of b^s < for the considered solutions inside the 
disk. The value of b^s can be determined by matching the 
calculated fields b r $ and b<ps onto an external field and 
flow as discussed in LRBK and here in §3. 

LRBK showed that there is no jump in b r across the con- 
ducting surface layer. This implies that du r /d(\Q s = 0, 
which represents a first boundary condition of the solu- 
tion of equation (8). The equations inter- reltaing b r , b$, 
u r , and Su^ are summarized in the Appendix. A second 
boundary condition, u r |^ s = 0, follows from equation 12 
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of LRBK evaluated just outside of the conducting layer. 
A third boundary condition, 9b ( f,/dC\^ s _ = 0, is derived in 
LRBK. 




\A_ 1.6 
u r /u 



Fig. 1. — Radial and toroidal field components normalized to 
B z (r, 0) at the disk's surface as a function of the average accretion 
speed u r (normalized by the viscous accretion speed Ho). For this 
plot p = 240 and V = 1. 



12 
10 

8 
6 
4 
2 







\ -n U r 1 

; \ u o i 






\ lml 


\v =1 




0.9 







200 220 240 260 280 300 320 

P 

Fig. 2. — Summary of solutions to equation (8) which satisfy 
the boundary conditions of §2.1. u r is the average accretion speed 
and uo is the accretion speed in the absence of an ordered magnetic 
field. The curves are smooth lines connecting the marked points. 
The physically relevant solutions have 1Z < 0[2\/27r/(a/3iio)] which 
is about 2.5 for a = 0.1, /3 = 240, and u = 0.15. 



3. SPECIFIC SOLUTIONS 

Here, to simplify the analysis we consider the limit 
where 7 — > 00 in equation (7) and S — > in equation 
(18) of LRBK. Then, (s — > Cm and both p and g are 
unit step functions going to zero at ( m — \J~2. Also, 
u r = (u r ) and X = Thus the above physical con- 

dition u r > 3efc„ = uo implies that b r s > uqCsV from 
equation (13) of LRBK. We assume k p = — d\np/d\nr 
and k v = 1. 

The solutions to equation (8) are u r oc exp(ifcj£) (with 
j = 1, 2, 3), where 



a 4 i3 2 (k]) 3 +2Va 2 /3{kf) 2 + {a 2 (3 2 +V 2 )k]-3l3V 2 = , (9) 



is a cubic in k 2 . The discriminant of the cubic is negative 
so that there is one real root, k\ , and a complex conjugate 
pair of roots, k\ 3 . Because u r is an even function of ( we 
can write 

u r = ai cos(fcoC) + CL2 cos(fc r C) cosh(fciC) 

+ a 3 sin(fc r C) sinh(fciC) , (10) 

where k = y/k^, k r = Re(-\/fcf), and h — Im(y / fcf). 

We consider a thin disk, e = h/r = 0.05, and a viscosity 
parameter a — 0.1. Figure 1 shows the dependences of the 
surface field components on the average accretions speed 
for j3 = 240 and V = 1. The b^s dependence is given by 
equation (6) and is independent of V while the b r s is given 
by b r s = V(su r from LRBK (equation 13). 



There are three boundary conditions at the surface of 
the disk. Additionally, there is the ratio 1Z of the accretion 
power going into magnetic disk winds to the viscous power 
dissipation. The solution for u r (C) has three unknown co- 
efficients, ax, 02- & a.3 which are dependent on a, "P, and 
p. If a and V are fixed, then there is a definite relation 
between 1Z and /3. Figure 2 shows this relation determined 
numerically. 




0.0 0.25 0.5 „ 0.75 1.0 1.25 1.4 



Fig. 3. — Radial inflow speed u r = —v r (normalized to ac s o) as 
a function of f = z/h for /3 = 240 and V = 1. The surface of the 
disk is at £ = \/2. 

Figures 3 shows the vertical profile of the accretion flow 
for /3 = 240. The profiles of Su^ = (v^ — vk)/vk, b r , b^ 
follow from the equations in the Appendix. 

The value of b c f,s+ < or TJ r /u is not fixed by the solu- 
tion for the field and flow inside the disk. Its value can be 
determined by matching the calculated surface fields b r s 
and b^s+ onto an external magnetic wind or jet solution. 
Stability of the wind or jet solution to current driven kink- 
ing is predicted to limit the ratio of the toroidal to axial 
magnetic field components at the disk's surface |6</,s+| to 
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values £ C(2tt) (Hsu & Bellan 2002; Nakamura, Li, & Li 

2007). Recall that TZ = u r /u - 1 = 2|6 0S+ |/ (aftV2u ) is 
the ratio of the accretion power going into the disk wind 
to the viscous dissipation in the disk. For the mentioned 
upper limit on |&0s + |, Tl £ O[2v^7r/(ay8u )] « 2.5. This 
implies a critical value of ft for a stationary disk solution 
(for V = 1), ft c w 1.2(ae)- 1 . For a = 0.1 and e = 0.05, 
& « 240. 

For ft ^ ft c , the disk and large-scale magnetic field is 
not be in a stationary state. Equation (6) for the av- 
erage accretion speed still applies and can be written as 
u r = uq + ub where ug is the magnetic contribution to the 
accretion due to the outflow of angular momentum from 
the disk's surfaces. In general, ub is an increasing func- 
tion of the magnetic field strength if TZ ~ 1 (Lovelace et 
al. 1994). For this reason for f3 < f3 c the ordered poloidal 
magnetic field threading the disk will be advected inward 
in the disk while for f3 > f3 c the field will advect outward. 
Thus, ft c acts as a threshold value for the buildup of a 
large scale field in the inner regions of an accretion disk. 

Another possibility is that the disk winds are negligi- 
ble so that TZ <C 1. Figure 2 suggests that in this limit 
we have a stable stationary solutions with ft values much 
larger than the mentioned critical value. 

4. CONCLUSIONS 

For a specific disk model with turbulent Prandtl num- 
ber V = 1 and TZ ~ l,we show that there is a critical value 
ft c w l.2r/(ah) for a stationary disk threaded by a large- 
scale magnetic field. For a = 0.1 and e = 0.05, ft c w 240. 
For ft < ft c , we argue that that the large-scale field will be 
advected inward and lead to the build up of the ordered 
field in the central region of the disk. For ft > ft c , the or- 
dered field will advect outward. Alternatively, if the disk 



wind is negligible so that TZ <C 1 , stationary disk equilibria 
with a large scale field are possible for ft much larger than 
the mentioned ft c . 

The 3D MHD simulations of Beckwith, Hawley, & Kro- 
lik (2009) show the inward advection of a significant frac- 
tion of the initial unipolar vertical magnetic flux thread- 
ing a MRI unstable plasma torus around a black hole. 
The initial average ft in the torus is 100 which is sig- 
nificantly less than the critical value estimated here. A 
broader range of simulations for much larger initial val- 
ues of ft would be needed to determine if the initial flux 
diffuses outward. Axisymmetric MHD simulations of Ro- 
manova et al. (2011) show episodes of enhanced accretion 
on the disk surfaces for the case of an MRI unstable disk 
around a rotating magnetized stars. However, the pres- 
ence of the stellar magnetic field complicates the analysis 
of flux transport in the disk. 

We have shown that the presence of a large scale poloidal 
magnetic field in an accretion disk around a black hole 
not only plays a decisive role in the jet formation, particle 
acceleration, and observed hard energy radiation (Bosch- 
Ramon 2011), but also creates the possibility for a new 
type of non-stationary behavior in these objects, which 
could be related to appearance of different spectral and lu- 
minosity states and transitions between them as observed 
in microquasars (Malzac and Belmont 2009). 
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APPENDIX 
SUBSIDIARY EQUATIONS 

Here, we summarize the relations between the field and velocity components. From the radial force balance, 

db r ftp h , 2 2\ , 2nd (_ du r \ 

where k p = —dlnp/d\nr is assumed positive and of order unity. The 0— component of equation (2) gives 

<% aftp , > aft d / duA 

= — ~Ur) ~ — ^ (flff-^ j , (A2) 

where k v = dln(pc^ r 2 /h)/d\n(r) > is also assumed to be of order unity. 
The toroidal component of Ohm's law gives 

— = -u r . (A3) 



The other components of Ohm's law give 



ditcf, 3e ae d ( db<p 



Combining equations (Al) and (A3) gives 
For thin disks, e <C 1, and /3 > 1, we have = 1 + Su^ with (du<p) 2 <C 1. Conseque 



. k v e 2 eVu r o?e d / du r 



to a good approximation. 



